First principles study on the spin transfer torques 
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An efficient first principles method was developed to calculate spin transfer torques in layered 
system with noncoUinear magnetization. The complete scattering wave function is determined by 
matching the wave function in the scattering region with the Bloch states in the leads. The spin 
transfer torques are obtained with aid of the scattering wave function. We applied our method 
to the ferromagnetic spin valve and found that the material (Co, Ni and Ni8oFe2o) dependence of 
the spin transfer torques could be well understood by the Fermi surface. Ni has much longer spin 
injection penetration length than Co. Interfacial disorder is also considered. It is found that the 
spin transfer torques could be enhanced by the interfacial disorder in some system. 

PACS numbers: 72.25.Ba, 85.75.-d, 72.10.Bs 
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I. INTRODUCTION 

Spin angular momentum can be transferred by the 
flowing electrons from one ferromagnetic (FM) material 
to another FM material, which is so-called spin transfer 
torques (STT) introduced by Slonczewskii and Berger— . 
Those two seminal studies have shown that the dynamics 
of magnetization in FM material could be dominated by 
the spin torques carried by electric current. The exci- 
tation of coherent precession of magnetization and spin 
wave were predicted. The STT was soon identified in the 
experiments^ by clear observation of the magnetization 
switching in FM spin valve, which excites great interests 
in experiment and theory^'^'^'^'^ i ' ^ ^ . 

The theories^'^'^ combining the quantum treatment 
of the interface scattering and the Boltzmann-like treat- 
ment of the bulk scattering work reasonable well with the 
experiments of metallic system. However, recent experi- 
ments on the tunnelling system^'^ and magnetic domain 
wall^ call for a full quantum treatment of the whole 
system. Edwards et.al.,^^ obtained the torques of spin 
valve in the empirical tight-binding frame and Haney 
et.al.f^ calculated the torques in the similar structure 
with nonequilibrium Green's function (NEGF) based on 
LCAO basis. 

Both semiclassical and quantum mechanical study 
show that the STT is most significant near the 
nonmagnet(NM)|FM interfaces in the spin valve. Up to 
now, only a few studies have addressed the material de- 
pendence of spin torque, which could be an important 
issue as the spin dependent transport is greatly affected 
by the electronic structure in FMi^ii^. Furthermore, pre- 
vious studies focused on ideal structure without consid- 
ering the disorder at the FM|NM interface, which should 
exist in the realistic spin valvei^. 

The main aim of this paper is to formulate a method 
to calculate STT of a noncoUinear magnetized system 
within the first principles frame. Differing from the pre- 
vious Green function based workii, we obtained the com- 
plete scattering wave functions of the whole systemic. 
The STT— is formulated in the tight-binding represen- 



tation. Large system such as domain wall can be well 
treated in this framework^-. We apply our formulism to 
the Co|Cu|FM|Cu spin valve system with impurity scat- 
tering at the FM|NM interface. Our study shows that 
the STT can penetrate deep into the ferromagnetic ma- 
terials for Ni, which is quite different from Co. It is also 
found that average torques are enhanced in the presence 
of interfacial disorder. 

This paper is organized as following. In Sec. II, we 
present the details of the formalism for constructing the 
eigenmodes of the lead and computing the STT in spin 
valve. Note that not only the transmission and reffection 
coefficient are obtained but also the wave function in the 
scattering regime is obtained explicitly. In Section III, 
the method is used to calculate the conductance and STT 
in the systems of Co|Cu|FM|Cu(lll), with FM is Co, 
Ni and Ni8oFe2o, respectively. The effect of interfacial 
disorder is discussed. In Sec. IV, we summarize our 
results. 



II. THEORETICAL MODEL 

Let us focus on the spin transport and STT in the 
layered systems sketched in FigH] The scattering region 
S, which is denoted by the layer index 1 < / < iV, is 
sandwiched by left(L) and right(R) leads. For this de- 
vice, there exists perfect lattice periodicity in the X-Z 
plane. Particle current flows along Y axis. In scatter- 
ing region no periodicity is assumed along current direc- 
tion. Here the atomic potentials were determined by the 
tight-binding linearized muffin-tin-orbital (TB-LMTO) 
surface Green's function (SGF) methodic. When com- 
bined with the coherent potential approximation (CPA), 
this method allows the electronic structure, charge, and 
spin densities of layered materials with substitutional 
disorder to be calculated self-consistently with high ef- 
ficiency. To model the noncoUinear system in the spin 
valve, the rigid potential approximation is used. In this 
approximation, we rotate the potential of fixed magnet 
in spin space to construct the relative angle between the 
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FIG. 1: (color online) Sketch of the configuration used for 
current-induced switching. A scattering region is sandwiched 
by left-(L) and right-hand(R)leads which have translational 
symmetry and are partitioned into principle layers perpendic- 
ular to the transport direction. The scattering region contains 
principle layers but the structure and chemical composition 
are in principle arbitrary. The switching layer FM can be Co, 
Ni, NigoFeao. 



the a representatioii2Si21 in the global coordinate system 
can be written as 



H 



xC/R,fe^y f^R']- (2) 



where Cr^^ and Ar^ are 2x2 potential parameter ma- 
trices expanded in spin space and diagonal in the local 
coordinate system. The unitary rotation matrix at site 
R can be defined by 



C^R (6'r, <^r) 



cos ■ 



, (3) 



polarization directions of fixed magnet and free magnet, 
which is a good approximation as the two magnets are 
spaced far enough by a Cu layer. 

Following previous work^^, we describe the theoretical 
frame developed with wave-function matching (WFM) 
based on TB-LMTO basis for studying the STT. In Sec. 
II A, we review the Hamiltonian and KKR equation for a 
device with noncoUinear magnetization. The equation of 
motion (EOM) for layered system is extracted from KKR 
equation. In the Sec. II B, the boundary conditions of the 
EOM are formulated in terms of the Bloch states in the 
leads. In the Sec. II C, by solving the EOM in the scat- 
tering region with embedding potentials of the two leads, 
we obtain the complete scattering wave function of the 
scattering region. In the Sec. II D and E, the particle 
current and spin current are formulated with those ob- 
tained scattering wave function expanded in TB-LMTO 
basis. 



A. Hamiltonian and KKR equation 

For layered systems, atoms can always be grouped into 
principle layers defined as to be so thick that the interac- 
tions between layers / and J ± 2 are negligible as shown 
in Figdl 

The EOM for /th principal layer can be written as 

H/,/_ia/_i + (H - E)^j ai + H/j+ia/+i = 0, (1) 

where E is always set to the Fermi energy Ep for the 
transport problem. Here, a/ is the a vector describing 
the amplitudes of the /th layer in terms of the localized 
orbital basis {RLQ, where R is the site index and L can 
be defined hy L = (l,m). I and m are the azimuthal 
and magnetic quantum numbers respectively. C =T (i) 
denotes that the basis is eigenstate in spin space, which 
is parallel (antiparallel) to spin quantization axis. 

To the first order approximation of the full LMTO 
Hamiltonian, a short-range TB-LMTO Hamiltonian in 



where 0r, (^r are the azimuth angles of the local quanti- 
zation axis. Screened structure constants Sr^ r,^, con- 
tain all information about the structure, which are block 
diagonal in the spin space. 



RL.R'L' 



*RL.R'L' 






*RL,R'L' 



(4) 



Note that Sr^ r,^, is spin independent. The Hamilto- 
nian of Eq. ([2|) yields eigenvalues corrected to first order 
in {E — Ep) and is exact when we set E — Ep. 

For a noncoUinear magnetized system, the "tail- 
cancellation" condition yields the KKR cquation^i, 

R'.L' 

where Cr^ — (Crl|,Crlx) has the relation to the 
wave amplitude aR^, of L orbital at site R as Cr^ — 

C/r ( Ar2^ 1 ?7^aRi . Pr^ {E) is the screened potential 

function matrix and contains all information about the 
atomic species at site R for calculating the electronic 
structure. It is diagonal in the local coordinate system. 



^RL (E) 



p'kl 

PrI 



(6) 



where Pr^ 



-",T(i) 



{e - Cr^ (Ar^^^' j and E is set 
to Ep for the transport problem we considered. 

As there exists two-dimensional translational symme- 
try in the lateral plane, the states along the transport 
direction can be characterized by a lateral wave vector 
k|| in the corresponding 2-dimensional Brillouin zone (2D 
BZ). The screened KKR equation in the mixed represen- 
tation of k|| can be expressed in terms of principal layers 



3 



S^y.Ci^i (k||) + (UiPij (Ef) uj 



Cr kii 



iC/+i(k||) =0, 



(7) 



where C/ (k||) is the wave vector describing the wave 
function amphtudes of the /th principal layer consisting 
of h atom sites and has the dimension of 2 (^max + 1) h ^ 

— ki 

2M. Pi J is 2M x 2M diagonal matrix, and Sj j is also 
2M X 2M matrices with its sub matrix defined 

by 



5, 



RL,R'L' 



Til 



exp [*k||.T||]5^ 



rl,(r'+T||)l' R'+T|| e /' ^ ' 



where / and /' are layer index and T|| is 2-dimensional 
translational vector in the plane of principal layer. 

Note that Eq.© is the EOM by analogy with Eq.dl]). 
We will solve it for a given energy Ep oi electrons to 
obtain the wave function of the scattering state. The ref- 
erence to k|| and Ep in the formulism will be suppressed 
in the following two parts Sec. II B and Sec. II C. 



B. Eigenmodes of the leads 



W(±) = C/oW(±) 

= c/o[wI(±). 



(±) 



w{ (±),... ,wi,(±)]. 
Following Ando2Si, define the Bloch factor as 



F{±)=W (±) A (±) (±) . 



where 
with 



A (±) is 
the diagonal 



(10) 



(11) 



a diagonal matrix 
elements given by 



A| (±),--- ,Ai, (±),A{ (±),... ,Ai,(±) 



local 



we have the 



In 

relation 



coordinate system 
C/(±) = F'^"'^(±)Cj(±>iI.It is easy to proof 
that the Bloch factor defined above satisfies the Bloch 
relation in global coordinate system 



C, (±) = F(±/~^C^(±). 



(12) 



Bloch factors matrix F (±) relates the wave amplitude in 
the /th layer to that in the Jth layer for a state in the 
lead. 



For the scattering problem, far enough away from the 
scattering region the wave function can be expressed rig- 
orously with asymptotic forms in terms of reflection and 
transmission coefficients and Bloch states in the leads. 
As the wave function should satisfy Bloch's theorem in a 
periodic potential, we set C„=A"Co. In local coordinate 
system, the EOM in lead becomes 



0,0) 



— 5'q ISi^o 




_Co 



_Co 



(9) 



Details for solving Bloch states Co can be found in 
Ref. 17]. To overcome the numerical difficult of the spin 
degeneracy in NM lead and reduce the calculation efforts, 
we solve the EOM in the leads in local coordinate sys- 
tem for each spin separately. In global coordinate system, 
Bloch states can be obtained after an unitary transforma- 
tion. For the amplitude of 0th layer, we have Co = UqCq. 

The propagating states and evanescent states can be 
identified and sorted into right-going(-l-) or left-going 
(— ). Letting v/J}^^ (±) denotes the solutions of C cor- 
responding to eigenvalue (±), where t(i) denotes the 
eigenstate parallel (antiparallel) to the local spin quanti- 
zation direction. Constructing the matrix W (±) as 



C. Scattering problem 

The equations of motion with open boundary con- 
ditions for a device usually contain infinite number of 
equations. By incorporating the boundary conditions 
in the leads, the scattering problem can be reduced to 
a set of coupled linear equations with finite number of 
equationsii. 

For an electron coming from the left lead, Eq.© for 
7 = can be rewritten as 



UoPo,oUq — S'o,oj Co — 5oaCi 
= 5o,-i[/^^-i(+)-/^^-^(-)]Co(+), (13) 

where L denotes the left lead and with 50,0 = <S'o,o + 
5o,-iF^^ (-). The Sq.-iFI^ (-) is the embedding po- 
tential for the left lead. 

In the right lead, only right-going waves exist in the 
(N + l)th layer. The EOM for / = A^ 1 is 

Vn+\Pn+\,N+\U^1<1j^\ ~ '5'Ar+i,jV+l I C^v+i — 5jv+i,ArCjv — 0, 



(14) 

where ^at+i^at+i = 5jv+i,Ar+i + S'Ar+i,jv+2^'i?, (+) and 
S-N+i.NJri^R (+) is the embedding potential for the right 
lead. 
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Making use of the lead boundary conditions for 0th 
and {N + 1) layer, the scattering wave function can be 



found as 
/ Co 
Ci 

C2 



Cat 



(upW- 



fS,,^, [F^' {+) - F^' H] Co (+)\ 



V 



, (15) 



where S is of block tridiagonal matrix containing Sjj 

except the ^o.o and Sn+i.n+i are defined as above. The 
spin polarization direction at different sites can be in- 
corporated by the unitary rotation U at corresponding 
site. 

To obtain the scattering state, we need to specify an in- 
coming state Co at the right side of Eq.ljTS]). This can 
be achieved by introducing the right going eigenmodcs of 
left lead as the incoming states by setting Cq to be 
w\ {+), where w\ {+) should be renormalized so as to 
carry an unit flux. Each corresponds to a scat- 

tering state in device. 

The amplitude of layers from to TV -|- 1 solved 
from Eg. p^ serves for computing the particle current 
and spin current. Also, the scattering matrix can be 
obtainedii. 



D. Particle Current 

Let us consider the particle current operator of a quasi 
one-dimensional TB model for a special k|| vector at E = 
Ep. The MTO-basis functions |RLC''") are obtained 
from the Bloch sum of the particle waves: 



|RiC''") = 



|R + T||,LC") 



(16) 



So the density operator at R site in the mixed represen- 
tation for a special k|| vectors can be defined by 



.k|| 
Pr 



^|RiC''")(R^C' 



(17) 



Neglecting the electron motion inside the atomic cells, 
the velocity operators can be expressed by the intersite 
hopping^'^ and will give the total current for subspace. 
The velocity (current) operator can be defined by 



1 

ih 



X,H 



(18) 



where X is the coordinate operator, which can be rep- 
resented in TB model by a diagonal matrix ^■r.l.r'L' = 

With aid of Eq.fTSl). the current operator Jr|/j^ from 
R'th to Rth site (R ^ R') can be written as 



<R'R I 



LL' 



r-Vk|| 

RL,R'L' 



h.C. 



(19) 



where H^'i^,R,i,<;, 
and H^'^ j^,^, is the Hamiltonian matrix in spin space, 
which has relation with Eq.([21) as 



Hi^'l^'L' = E^-P • T||] H^,,(^,+T„)^,. (20) 

Til 



The expectation value of operator A is <^A 
A The particle current can be expressed as 

Jr'r (k||)) = E^t^Ri (k|l) ^nL,ii'L'^R'L' (k||)-/i.c.], 

LL' 

^ (21) 

where where brl (ky) = (hrlt (k||) ,aRLj, (k||)) , and 
aRLC (k||) = (RiC"*"" I*)- aRL (k||) has the relation with 
Crl (k||) as follow 

aRL (k,|) = C/r (Ar^)"" uICrl (k||) . (22) 

The Crl (k||) can be obtained by Eq. p5|) for a given 
fc|| . Within the MTO formulism, the current can also be 
expressed with structure constants matrix as in RefflTl 

Jr'R (k||)) = Eit^RL (k||) 5^^.R'L'Cr,z.' (k||)-/..C.] 



LL' 



(23) 

The continuity equation of particle current at R site 
in the /th principle layer reads 



R'e/-i,/ 



E 



R'e/,/+i 



Jr.r' (k 



dt 



(24) 



where the first term at the left side of Eq. ((24|) is the 
incoming current to the R site and the second term is 
outgoing current from this site. 

As shown in Fig. ([21), the current is assumed to flow 
from / — 1th layer to / -|- 1th layer. Considering the R 
site, the incoming current is composed of current from 
the sites in the / — 1th layer and the sites ahead R site 
(relative to transport direction) in /th layer. If there ex- 
ist other atoms in the same plane of R(see Fig. ([2])), the 
current from those atoms also are considered as the com- 
ponent of incoming current to R site. The outgoing cur- 
rent is composed of current to the sites in / -I- 1th layer 
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1-1 



I+l 



FIG. 2: Illustration of incoming current and outcoming cur- 
rent for the Rth site. Assuming the particle current comes 
from / — 1th layer to 7 + 1th layer. Arrow lines denote the 
current related to Rth site and dot lines denote the coupling 
between sites irrelevant to Rth site. 



and those sites behinds R site in Ith layer. Note that 
treating the current between atoms in the same plane as 
the incoming current or outgoing current will not result 
in any physical consequence. Careful check has been car- 
ried out that the particle current conservation law can be 
satisfied atom by atom and layer by layer. For the scat- 
tering states we calculated, the right side of the Eq. ((Ml ) 
is zero. 

In the linear response regime, the particle current un- 
der a small bias Vb at zero temperature can be expressed 



'RR' 



e 

hNi 



-E 



'RR' 



(25) 



where the bias is given by the difference between the 
electrochemical potentials of the two leads as eVf, = /i^ — 
fiR, and A^ii is the number of kn points in 2D BZ. 



E. Spin current and spin torque 

The spin current is defined similar to the particle cur- 
rent in Section II. D. Considering a quasi one-dimensional 
TB mode for a special ky vector, the spin density opera- 
tor at site R can be defined as 



'R 



(26) 



where is 2 x 2 Pauli spin matrix. The spin current 
operator generally can be defined as 



o-(g)V + 



(27) 



note that is a tensor. For spin current between Rth 
and R'th site (R ^ R'), we could project J along the 
direction vector xr_r' in real space as j • xr_r' . Then 
the spin current operator Jr'^r (ky) from R'th to Rth 
site (R 7^ R') can be written as 



Jr',R (k||) = Y.^^^^'^L.^'L' + H^L,R'L'^ - h.C.]. 

(28) 



LL' 



where Jr'.r (ky) is a vector only in spin space. 
For a specific state j^f), the expectation value is 



Jr' 



R I 



LL' 



T- 



^RL 



(k„)aH 



RL.R'L 



R'L' 



r 



^RL 



k,i H: 



RL.R'L'C^^R'i' 



kii -/i.. 



(29) 



The STT (T^ (k||)^ can be defined as the difference 

of the incoming spin current and outgoing spin current 
of R site in the Ith principal layer: 



Tr k|| 



E (^R'3(k|l))- E (^R,R'(k 



R'e/-i,/ 



R'e/,/-i-i 



(30) 



where the superscript s is used to denote the incoming 
state is parallel or antiparallel to the local spin quantiza- 
tion axis of injection lead, which is very helpful, e.g. in 
ferromagnet we can distinguish the contribution to the 
total torques from the majority spin or minority spin. 
Such definition consists with those in Ref.fll], where an- 



alytic analysis shows that for STT defined in this way 
equals to the exchange torques acted on the injected spin 
defined in Eq.f^H) with only a sign difference. After sum- 
mation over 2D BZ, spin torque acted on R th atom can 
be expressed as 



-R 



(31) 



where the bias is given by the difference between the 
electrochemical potentials of the two leads as eVf, = Ml — 
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FIG. 3: (color online) (a)The total conductance of 
Co(6')|Cu(9ML)|Co(15ML)|Cu versus polarization direction 6 
of fixed magnet Co. (b) The angular dependence of total spin 
torques on free magnet Co, where the electron current flow 
from the fixed magnet to the free magnet. 



III. SPIN TRANSFER TORQUES IN 
CO/CU/FM/CU (111) SPIN VALVES 

A. Ordered interfaces 

A spin valve of Co|Cu|FM|Cu as shown in Figgis used 
as an example to illustrate our method. The left lead 
consists of semi-infinite Co with the polarization direc- 
tion 9 (see FigH]). Cu spacer of 9 monolayer (ML) is 
located between fixed magnet Co and free magnet FM. 
The free magnet contains d ML, which could be Co, Ni, 
or Ni8oFe2o in this study. The lattice constants is as- 
sumed to be uniform in the whole spin valve, that is, 
acu = ofm ~ 3. 54 A and the transport is along /cc [111]. 
With spd -basis, exchang-correlation potential is calcu- 
lated and parameterized by Vosko-Wilk-Nusair— . Our 
calculation gives the magnetic moments as 1.64/^s/Co 
atom, 2.60/ZB/Fe atom and 0.62/iB/Ni atom. For the 
calculation of transport, total 90000 k|| points in 2D BZ 
are summed. 

Firstly, we present the angular dependence of total 
conductance G{9) of the spin valve with the free mag- 
net to be 15ML Co in Figl3] (a). The monotonic de- 
crease with increase of 9 is consistent with the previous 
ah. initio results--. Giant magnetoresistance (GMR) can 
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Layer index 

FIG. 4: (color online) The layer dependence of STT on in- 
terfacial unit cell in the spin valve Co{6 = 90)|Cu|FM|Cu , 
where the free magnet FM are Co, Ni, Ni8aFe2o respectively. 



be defined as GMR \ 



G(0°)-G(180°) 



100%, which is 24% 



G(180°) 

in this case. With electron current flowing from the fixed 
magnet to the free magnet, Fig[3] (b) gives the angular 
dependence of total spin torques on the free magnet Co, 
which restore the line shape of spin torques obtained in 
previous workii. With the drive of the in-plane torque, 
magnetization of free layer is going to parallel to that of 
fixed layer. Due to the breakdown of time reverse symme- 
try for spin current, if the direction of electron current is 
reversed, the in-plane torques is going to drive free layer 
to antiparallel to fixed layer. Such phenomenon is exactly 
the current induced switching of magnetization observed 
in spin valve. 

Layer resolved Spin Torque: The layer resolved STT 
contains the information about whether the spin angle 
moment is absorbed near the interface or not. FigUgives 
the comparison of the layer dependence of STT in the 
spin valves with three different free magnet. Here the 
polarization of the fixed magnet Co is set to 6* = 90° 
without lost of generality. In our frame, Ta; corresponds 
to the in-plane torques and Ty corresponds to the out- 
of-plane torque. The decay and oscillation of the STT 
are greatly different among those materials we studied. 
When free magnet is Co as shown in Fig|4{a), our re- 
sult almost reproduces the previous resul t^°i^^ . The fast 
decay of the STT indicated the surface atoms absorbed 
most of the spin angle moment as the current passes by. 

However, when Ni serves as the free magnet as shown 
in Fi^Ub), the maximum torques is not on the surface 
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FIG. 5; (color online) The layer dependence of STT on intcr- 
facial unit cell for the spin injection setup of single interface, 
(a) Cu|Co (b)Cu|Ni. 



FIG. 6: (color online) (a) and (b) The layer dependence of 
STT on interfacial unit cell when the spin injected through a 
single Cu|Co interface for different ky points in 2D BZ. 



atom and the decay is very slow with much longer oscil- 
lation. This observation is quite different with our pre- 
vious knowledgoi^. The similar behavior is also found in 
Fi^Jc) as Ni8oFe2o is free magnet, the oscillation looks 
like that in Ni, but decaying faster. Due to the lack of 
obvious oscillation, the total in-plane torques on NigoFe2o 
(7.7 X lO-^eVh) is greater than that on Co (5.0 x IQ-^eVb) 
and that on Ni (6.3 x IQ-^eVt). 

The layer resolved STT shown in Fig|4] could be af- 
fected by the multiple scattering between the two inter- 
faces with Cu. To remove multiple scattering effect on 
the torque, we perform the calculation for single inter- 
faces of Cu(90°)|Co and Cu(90°)|Ni, with 100% polar- 
ized electrons injected from Cu side. Here Cu(90°) indi- 
cates the polarization direction of the injected electrons. 
The results are shown in FigO For both interfaces, the 
maximum torques exists on the surface atoms and the 
oscillation spectra in ferromagnet become smoother and 
clear. Still the oscillation behavior strongly depends on 
the materials. The STT exists only near the Cu|Co inter- 
face, while the STT penetrates deep into the Ni for Cu|Ni 
interface. Due to the long penetration length, the multi- 
ple scattering between two Cu|Ni interface does appear 
in FigH^b). 

Simple Model for Spin Torques in FM:Fot the layered 
system such as spin valve, the incoming state of the in- 
jection lead can be labelled by k|| in 2D BZ. Generally, 
these states will be coupled to the propagating states 
and evanescent states of another side of the injection in- 
terface. The STT can be expressed as 

r OC ^ C^^e*[('=i-'=i)^+'^-^l + 3decaj;(x), (32) 



where first term denotes the contribution from the prop- 
agating states^ and fc^ — kl gives the spatial precession 
frequency Ak^^,. The contribution from the propagating 
states should oscillate as function of position and will 
not decay as shown in FiglSJ^b). However, the frequency 
could be quite different as k|| runs over the 2DBZ, so the 
final contribution will decay after summation over 2DBZ. 
The second term of Eq. ([5^. Qdecay{x), is the contribu- 
tion from the evanescent states. As we have known that 
no particle current can be carried by evanescent state, 
however, such states do give effect on the spin current 
and also on the STT. Evanescent states do contribute 
to spin torques and should responds for the initio decay 
of the STT in the system as Co|Cu|Co|Cu, as shown in 
Figl6^a) where the evanescent state dominates. 

Two reasons could account for the decay of the STT 
away from the interface, (i) Vanishing of the evanescent 
states' contribution. For Cu|Co, our calculation shows 
that this part of contribution is about 10% of the total 
torques on the first layer close to injection interface, (ii) 
Cancellation effect among different k|| in 2D BZ^. 

The materials dependency of the STT could be under- 
stood based on the Fermi surface of Ni and Co. The wave 
vector fc^(^) can be found by the projection of minority 
spin (majority spin) Fermi surface of ferromagnet along 
the current direction, where fi (v) denotes the different 
sheets of Fermi surface for minority spin (majority spin). 
In Fig[7l the Fermi surface for Co and Ni viewed along the 
(111) direction for majority and minority spins is shown. 
As the shape of Fermi surface for majority spin and mi- 
nority spin in Co are greatly different to each other, the 
precession frequencies Ak^^, of injected spin are varied 
rapidly as k|| running over the 2DBZ. After summation 
of 2D BZ, the strong cancellation is expected as shown 
in FiglU^a) and FigOJa). However, for Cu|Ni, due to 
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FIG. 7; (color online) The First row; Fermi surface projection 
of the Co bulk fee Brillouin zone on to a plane perpendicular 
the (111) direction. Left-hand panel is for majority electron 
and right-hand panel is for minority electron with band 3,4,5 
FS. The Second row is Ni bulk. 



the similar symmetry between the wave function of the 
sheet 11 = 6 oi minority spin and that of Cu,the electrons 
pass Cu|Ni interface mainly through this channel(/i = 6). 
While for majority spin, the Fermi surface contains only 
one sheet. The precession frequency is dominated by 
A/cge- Similar shape of sheet /i = 6 of minority spin 
and sheet ^ = 6 of majority spin will result in amount 
of propagating states with similar precession frequency. 
After summation of those states, the cancellation could 
be week and collective oscillation must be of long period. 

The above physical picture about spin torques in FM 
should be qualitatively applicable to the system with FM 
to be Ni8oFe2o. For Ni8oFe2o, the overall band structure 
resembles that of Ni, however, due to the scattering of 
Fe impurity atoms, the fine structure at Fermi surface 
could be much more complicated than that of Ni. The 
dispersion of precession frequency A/c^i, is large and the 
cancellation should be strong. As a result, the decay of 
spin torque is much faster than in the conventional FM, 
Ni and Co. Besides, spin-orbit couphng is not included 
in our calculation yet, which could introduce new mech- 
anism of decay in Ni8oFc2o- 



B. Interfacial disorder 

Interfacial disorder is likely to exist in the metallic sys- 
tem. Previous studiesi^ showed that the interfacial alloy 
could change the polarization of the interface resistance. 
How the interfacial disorder affects the STT is question 
we would like to answer in this section. The interfacial 
disorder is introduced by two layers of substitutional al- 
loy Cua;FMi_2; and Cui_a;FMj;at interface between Cu 
and FM. In present study, alloy is modelled by a 8 x 8 
lateral supercell, which was shown to be a good modelling 
of the interfacial alloyii. 

In Fig[Ha)&(b), the concentration x dependence of 
total conductance and total torque(in-plane and out-of- 
plane) on free magnet are given for the realistic spin valve 
of 00(6* = 90°)|Cu|Co(Ni)|Cu. The total conductance de- 
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FIG. 8: (color online) The concentration x dependence of the 
conductance (a) and total spin torques on the free magnet 
(b) of spin valve Co(6l = 90)|Cu|FM|Cu with interfacial alloy 
Cu^FMi-a;, where the solid symbol for FM to be Co and open 
symbol for FM to be Ni. In (b) , the black symbol for in-plane 
torques and red symbol for out-of-plane torque. 
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FIG. 9; (color online) The layer dependence of STT on inter- 
facial unit cell in the spin valve Co{6 = 90)|Cu|FM|Cu with 
interfacial disorder, (a) FM is Co, (b) FM is Ni. The zone la- 
belled by n is the layers with substitutional alloy CusoFMso. 
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FIG. 10: (color online) The thickness d dependence of the 
averaged spin torques on the atom in free magnet and the 
total conductance of the spin valve. 



creases with increase of the concentration which means 
the interfacial disorder will suppress electronic transport 
in this system. However, as shown in FiglHlJb), the total 
in-plane torques acted on free magnet increase when the 
disorder is enhanced in spite of the almost constant out- 
of-plane spin torque. The increase is around 50% for Co 
and 30% for Ni. 

In Figini the layer dependence of spin torques is shown. 
Comparing with FigHJa), it is found that the torques on 
atoms near the interface with Cu spacer have been en- 
hanced. This result means that the interfacial roughness 
will not kill the interfacial spin torques, on the contrary, 
the dirty interface may be helpful to enhance the torques 
transfer. 

For the spin valve with FM is Co, the free magnet 
thickness d dependence of the total conductance and STT 



are shown in Fig llOi where the spin torques is obtained 
by average of total torques on free magnet over all atoms 
in free magnet. Due to the quantum size effect, the con- 
ductance decays with a small oscillation and tend to be 
constant with increase of free magnet thickness. The in- 
plane torques dominates over the out-of-plane torques for 
all thickness. 



IV. SUMMARY 

Based on the first principles frame, a method was de- 
veloped to calculate the transport and spin torques of 
the layered system with noncoUinear magnetization in 
linear response regime. STT in the ferromagnetic (FM) 
spin valves are calculated. We found that the behavior 
of spin torques in the free layer are greatly dependent 
on the materials. The cancellation effect of the STT due 
to the different precessional frequency is sensitive to the 
band structure of material. The contribution of evanes- 
cent states to the STT is found to be nontrivial at the 
NM|FM interface. The effect due to interfacial disorder 
is also considered, it is found that average torques are 
enhanced in the presence of disorder. 
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